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Abstract 

The transition form factors of D s — > foil', D — > /o-fe' and -B u — > /o^f decays are 
calculated in 3— point QCD sum rule method, assuming that /q is a quark-antiquark 
state with a mixture of strange and light quarks. The branching ratios of these decays 
are calculated in terms of the mixing angle. 
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1 Introduction 



The inner structure of the scalar mesons in terms of quarks is still an open question in par- 
ticle physics and it is the subject of an intense and continuous theoretical and experimental 
investigations for establishing their their nature (for a review, see [1]). There are numerous 
scenarios for the classification of the scalar mesons. The established ++ mesons are divided 
into two groups: 1) Near and above 1 GeV, and 2) in the region 1.3 GeV -r 1.5 GeV. The 
first group scalar mesons form an SU(3) nonet, which contains two isosinglets, an isotriplet 
and two strange isodoublets. In the quark model, the flavor structure of these scalar mesons 
would be 

|<7 > = cos 9\nn > — sin 9\ss > , 
\fo > = cos 9\ss > + sin9\nn > , 

a° = —/= (uu — dd) , <2q = ud , = du , 

k + = su , R° = ds , k~ = us , (1) 

where \nn >= (uu + dd)/y/2, and 9 is the mixing angle. Here we take into account the fact 
that between isoscalars ss and uu + dd there is mixing, which follows from experiments. 
Indeed the observation 

r(j/v> - M ~ \v(j/^ - M) , 

indicates that the quark content of / (980) is not purely ss state, but should have non- 
strange parts too [2]. Secondly, if / (980) is purely ss state, then f — > mi should be 
OZI suppressed. But the decay width of /o(980) is dominated by /o — > hit which leads to 
the conclusion that in /o(980) there should be fin parts as well. Therefore fo should be a 
mixture of ss and fin, as is presented in Eq. (1). Analysis of the experimental data shows 
that the mixing angle 9 lies in the range 25° < 9 < 40° and 140° < 9 < 165° [3]. 

Although there is another scenario where mesons below or about 1 GeV is described as 
a four-quark state (see for example [4]), in this work we restrict ourselves to considering the 
qq description for / (980) meson, but taking into account the mixing between ss and nn. In 
the present work we study the semileptonic decays B + — > / (980)£ + z/, _Dj~ s — > / (980)£ + z/ 
decays in order to get information about the quark content of / (980). 

From theoretical point of view, investigation of the semileptonic decays is simpler com- 
pared to that of hadronic decays, because leptons do not participate in strong interactions. 
The experimental study of weak semileptonic decays of heavy flavored mesons is very im- 
portant for the more accurate determination of the Cabibbo-kobayashi-Maskawa (CKM) 
matrix elements, their leptonic decay constants, etc. 

The precise determination of the CKM matrix elements depends crucially on the possi- 
bility of controlling long distance interaction effects. So, in study of the exclusive semilep- 
tonic decays the main problem is calculation of the transition form factors, which involve 
the long distance QCD dynamics, belonging to the non-perturbative sector of QCD. For 
this reason, in calculation of the transition form factors some kind of non-perturbative 
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approach is needed. Among all non-perturbative approaches QCD sum rules method [5] 
is more powerful, since it is based on the first principles of QCD. About the most recent 
status of QCD sum rules, the interested readers are advised to consult [6]. 

Semileptonic decays D -> K°eP e [7], D + -> K(K°*)e + is e [8], D -> neu e [9], D -> pez> e 
[10], 5 — > D(D*)£u£ [11] and £> — > 0£z/£ [12] are all studied in the framework of 3-point 
QCD sum rules method. 

In this work we study the semileptonic B u — > / (980)£ + z/^ and D s u) ~^ /o(980)£ + z^ 
decays in the 3-point QCD sum rules method. The paper is organized as follows: In 
section 2, we derive the sum rules for the form factors, responsible for pseudoscalar to 
scalar meson transition. Section 3 is devoted to the numerical analysis of the transition 
form factors and discussion and contains our conclusions. 



2 Pseudoscalar-scalar meson transition form factors 
from QCD sum rules 

Pseudoscalar-scalar transition form factors are defined via the matrix element of the weak 
current sandwiched between initial and final meson states (S(p') 1^17^(1 — 75)52! P(p)), 
where 51 and q 2 are the relevant quarks, P and S are the pseudoscalar and scalar me- 
son states, respectively. It follows from parity conservation in strong interaction that only 
axial part of weak current gives non-zero contribution to this matrix element, and imposing 
Lorentz invariance, it can be written in terms of the form factors as follows: 

(Stf) |^(1 - 75)g 2 | Pip)) = -i[f+(p + P% + f-q»] , (2) 
where g M = p 1 - p 2 . 

For evaluation of these form factors in the QCD sum rule, we consider the following the 
3-point correlation function 

n>V 2 , q 2 ) = - I d 4 xd 4 ye^'y-^ (0 \t{ J s (y) J* (0) J P (x)}\ 0) , (3) 

where, J5 = 5252, J a = <?27^75<?i and Jp = 517552 are the interpolating currents of scalar 
and pseudoscalar mesons, and weak axial currents, respectively. It should be noted here 
that, 53 = u, 52 = u and 51 = b for the B u — ■> / (980) transition; and 53 = s(d), q 2 = s(d) 
and 51 = c for the D s ^) — > /o(980) transition, respectively. 

The decomposition of the correlation function (3) into the Lorentz structures, obviously, 
has the form 

n M = n + (p + P 'v + n-(p-p'V • (4) 

For the amplitudes 11+ and Il_, we have the following dispersion relation 

n±(pV 2 , Q 2 ) = / ^T^73^ + subtraction terms , (5) 

where p± is the corresponding spectral density and Q 2 = —q 2 > 0. According to QCD sum 
rules approach, the correlation function is calculated by the operator product expansion 
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(OPE) at large Euclidean momenta p 2 and p' 2 on one side, and on the other side it is 
calculated by inserting a complete set of intermediate states having the same quantum 
numbers with the currents Js and Jp. 

The phenomenological part of (3) is obtained by saturating correlator it with the lowest 
pseudoscalar (in our case B u , D s or D mesons) and scalar / (980) mesons, yielding 



(0\J s \S(p'))(S(p') #(0) P(p))(P(p)\J P (x)\0) 



(m 2 s — p' 2 ) (m 2 p — p 2 ) 
The matrix elements in Eq. (6) are defined as 

(0\J s \S(p')) = Xs, 
(P\J P \0) = -i 



+ excited states . 



(6) 



mpf f 



mi + m 2 



(7) 



where fg and fp are the leptonic decay constants of scalar and pseudoscalar mesons, and 
ms and mp are being their masses, respectively. Note that, leptonic decay constant fs in 
Eq. (7) is scale dependent for which we choose the scale to be \x = 1 GeV 2 , and 



mi 



m b for B u -> fo£v , 

foiv , D 



m 



for D s 



fo 



m 2 



m 



„ for B u 



foils, D — > foils , 

m s for D s -> /o^ , 
Using Eqs. (2), (4), (6) and (7), for the invariant structures we get 

fpm 2 P X s f± 



n 



± 



mi + m 2 (m| — p' 2 )(m 2 P — p 2 ) 



(8) 



From QCD side, the correlation function can be calculated with the help of the OPE 
at short distance, and in this work we will consider operators up to dimension six. The 
theoretical part of the correlator for the B s — > D so (2317)iis is calculated in [13], and in the 
present work, for the theoretical part of the corresponding sum rules, we will use the results 
of this work. 

For the spectral densities we have 



n+ = 



n_ = 



4A 1 /2( SjS ' ) Q2) 

Nr. 



4A 1 /2( S)S ' ) Q2) 

+ A' - A - 2m 1 m 2 



(A' + A) (1 + A + B) + (ml + 2 mi m 2 + Q 2 )(A + B) 
(A' + A + m\ + 2mim 2 + Q 2 + 2m l m 2 ") (A - B) 



(9) 



(10) 
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where N c = 3, A = s — m 2 , A' = s' — m 2 ,, and 



A = 
B = 



A( S ,s',Q 2 ) 
1 

A(*,s',Q 2 ) 



- (s + s' + Q 2 )A' + 2s' A 
-(s + s' + Q 2 )A' + 2sA' 



For the decays under consideration, m 2 is m u (m,d) or m s , and therefore, to take into account 
SU (3)-violating effects, here and in all following calculations we will retain terms that linear 
with 777.2, while neglecting the terms higher order in m 2 . 
For power corrections (PC)we get 



t 



m 1 — m 2 1 



rr 



+ ~;rn 2 {q 2 q 2 ) 



mf 



rr' 



1 

12 



m 2 q 2 q 2 ) 



?>m\(mi — m 2 ) 2(?77i — 2m 2 ) 2(2?77i — m 2 ) 

' /2 ' 



ry^t ryf 



+ 



mi (2m 2 + mim 2 + 2Q 2 ) — 2m 2 (m 2 + Q 2 



^* 2 ^* /2 
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+ -^Tra s (q 2 q 2 ) 2 



12777^(7771 — ^2) 8mi?77 2 (777 2 + Q 2 ) 56?77im 2 

— 1~ o .0 



ry^.ryf 



ry^ ry ^3 



/3 



Am\{2m\ + mim 2 + 2Q 2 ) — 8mim 2 (m 2 + Q 2 ) 87771(87771 — 7m 2 ) 



48 48 

+ ^7 + 



4(5m 2 - 20m 1 m 2 - 2Q 2 ) 



rr 
1 



/2 



^ 2^/2 



+ ^ra 2 m 2 (q 2 q 2 } 2 



m 2 (m 2 + Q 2 ) hm\ + AQ 2 Qmf 10m 2 



ry^ry ^2 



+ 



^ 2^/2 



+ 



+ 



(11) 



_ pr 1, .m 1 +m 2 1 / mi 

n_ = --(^2) ; + - A m x m 2 [q 2 q 2 ) - 

A rr 4 \ 



+ ^ m o^2> 



3m 2 (mi + 777 2 ) 2(mi + 3m 2 ) 2(3mi + m 2 ) 



rr 



/2 



+ 



mi (2m 2 + mim 2 + 2Q 2 ) + 2m 2 (m 2 + Q 2 ) 



ry 2y»/2 



12mf(mi + m 2 ) 8mim 2 (m 2 + Q 2 ) 56mim 2 



ry^ryt 



ry*2i ryt^t 



/3 



Am\{2m\ + mxm 2 + 2Q 2 ) + 8m 1 m 2 (ml + Q 2 ) 8mi(9mi + 7m 2 ) 



ry^t ryf^ 



ryt^S ryf 



28m 2 
+ -^-=r + 



ry^ryf 1 ^ tytyf^ /y*2i/yf 



+ gmom 2 (g 2 g 2 ) 2 



m 2 (m 2 + Q 2 ) 



qy^>ry! ! 2i 



ly^/yf^ 



6mi 4 
1 + 



24m 2 ' 



ry^ryf ryrylz\ ry^iyf ry^iy/ 



(12) 
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where r = p 2 - m\ and r' = p' 2 . Note that the D s -> f (980)l + is e and L> -> / (980)£+^ 
decays which are considered in [14] differ from our results in three aspects: 

• Our result on spectral density is two times smaller compared to that given in [14]. 
Since it is known that the main contribution to the sum rules comes from the spectral 
density, it is indispensable that our results on the form factors differ from those 
predicted in [14]. 

• In [14], part of those diagrams which are proportional to m s are not taken into account 
(in our case they correspond to the terms proportional to m2mo(<?2<?2))- 

• Sum rules for the form factor /_ are totally absent in [14], which could be essential 
for the B u — > f (980)ris T decay. 

Contribution of higher states in the physical part of the sum rules are taken into account 
with the help of the hadron-quark duality, i.e., corresponding spectral density for higher 
states is equal to the perturbative spectral density for s and s' Q starting from s > s Q and 
s' > Sq, where s and s' ar the continuum thresholds in the corresponding channels. 

Equating the two representations for the invariant structures IT-|-, and applying double 
Borel transformation on the variables p 2 and p' 2 (p 2 — > M 2 , p' 2 — > M' 2 ) in order to suppress 
the higher states and continuum contributions, we get the following sum rules for the form 
factors /+ and /_: 

= _ m 1 +m 2 l ml/M2em y M/2 ( r ( J Q2 )e -./M*-S,M» 

jpm A P \ s [J 
+ B M 2B M ,2U^ . (13) 

The double Borel transformation for the quantity \jr n r' m is defined as: 

r n r im Y[n) T(m) 

The integration region for the perturbative contribution is determined from the following 
inequalities: 

2ss' + (m 2 -s)(s + s' + Q 2 ) 
~ \V 2 (s,s',Q 2 )(m 2 -s) ~ { } 

In the calculation of the widths of the considered decays, it is necessary to know the q 2 
dependence of the form factors in the whole physical region m 2 < q 2 < q 2 



imax' 



3 Numerical analysis 

In this section we present our results for the form factors f+(q 2 ) and f-(q 2 ) for the decays 
under consideration. The main input parameters for the sum rules are the Borel parameters 
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M 2 and M' 2 , continuum thresholds So and s' . The values of other parameters needed are: 
m b = (4.7±0.1) GeV [6], m c = 1.4 GeF, m s = 0.15 GeV, (uu)]^ GeV = -(0.243) 3 GeV 3 , 
{ss) = 0.8 x (uu) [15]. The values of the leptonic decay constants of B u , D s and D 
mesons are determined from the analysis of the corresponding two-point correlators: Jb u = 
(0.14±0.01 GeV [16], f Da = (0.22±0.02 GeV [17] and f D = (0.17±0.02 GeV [6, 16, 17]. For 
the continuum thresholds we take the values Sq u = (33 ±2) GeV 2 , Sq s = (7.7 ± 1.1) GeV 2 , 
sff — (6 ± 0.2) GeV 2 and Sq = 1.6 GeV 2 which is determined from 2-point sum rules 
analysis [6, 14, 16, 18]. 

The Borel parameters M 2 and M' 2 are the auxiliary parameters and therefore the phys- 
ical quantities should be independent of them. For this reason we need to find the working 
regions of M 2 and M' 2 where form factors are practically independent of them. 

In obtaining the working regions of M 2 and M' 2 the following two conditions should be 
satisfied: 

• The continuum contribution should be small, and, 

• power corrections should be convergent. 

Our numerical analysis shows that, both conditions are satisfied in the region 10 GeV 2 < 
M 2 < 20 GeV 2 for B u -> f £i? e , 4 GeV 2 < M 2 < 8 GeV 2 for D S (D) -> f £u e , and 
1.2 GeV 2 < M' 2 < 2 GeV 2 for all channels. 

Varying the input parameters sq, s' , fo, fr> 3 , /# and f D in the respective regions as 
mentioned in the text, we get the following results for the form factors at q 2 = 

1.7(0.25 ±0.02) , 
-1.7(0.24 ±0.03) , 
1.7(0.32 ±0.03) , 

1.7(0.27 ± 0.02) , (16) 

The multiplying factor 1.7 corresponds to the case for As = 0.19 GeV 2 , and without this 
factor As = 0.35 GeV 2 [3]. For a comparison we present the results of the form factor f + 
for the B u — > / (980) transition coming from the covariant light front dynamics [19] and 
dispersion relation approach [20], as well as the result for the B u — > n transition [21]. 

f+ u ^ h (0) = 0.27 [19] , 
/+^ /o (0) = 0.09 [20] , 

fB»^*(0) = 0.25 [21] . (17) 

From a comparison of Eqs. (16) and (17) we see that, our prediction on / + for the B u — > 
/o(980) transitions quite close to the prediction of the light front dynamics and that of 
B u — > 7r transition when As = 0.35, and approximately three times larger compared to that 
of the dispersion relation approach. These close results of the form factor /+ for the B u — > n 
and B u — > f transitions could indicate of the that As should have the value As = 0.35, 
which is obtained in [3] by taking 0(a s ) corrections into account. 

Note that we present the form factor /_ only for the B u — > j^rv T decay, because this 
form factor can give considerable contribution to this decay. 



/+"(0) = 
/f.(0) = 

/?(0) = 
/?*(0) = 
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In estimating the width of P — > fo(980)£i?£ decay, we need to know the q 2 dependence of 
the form factors f+{q 2 ) and f~{q 2 ) in the whole kinematical region m 2 < q 2 < (mp — mf ) 2 . 
The q 2 dependence of the form factors can be calculated from QCD sum rules (see [8,9]). 
Unfortunately QCD sum rule cannot reliably predict q 2 dependence of the form factors 
in the full kinematical region. The QCD sum rules can reliably predict q 2 dependence 
of the form factors in the region approximately 1 GeV 2 below the perturbative cut. In 
order to extend the dependence of the form factors on q 2 to the full kinematical region, 
we look such a parametrization of the form factors where it coincides with the sum rules 
prediction of in the above-mentioned region. Our numerical calculations shows that the 
best parametrization of the form factors with respect to q 2 are as follows: 



/p(0) 



(18) 



1 — a P q + bpq 2 — Cpq 3 + dpq 4 ' 

where P = B u , D s , D and q = q 2 /m 2 P . The values of the parameters fp(0), ap, bp, cp and 
dp at As = 0.19 GeV 2 , are given in table I. 





/ + (o) 


/-(0) 


a 


b 


c 


d 


D s 


1.7 x 0.27 




0.87 


-0.17 


0.37 


1.46 


D 


1.7 x 0.32 




0.89 


-0.40 


0.18 


-1.00 


B u 


1.7 x 0.25 




0.48 


-0.30 


-0.47 


-0.99 


B u 




-1.7 x 0.24 


0.41 


-0.42 


-0.95 


-1.55 



Table 1: Form factors for the D s — > f £i/g, D — > f £vg and B u — > f £vi decays in a four- 
parameter fit. 

The dependence of the form factors / + and /_ (for B u — > f§Ti> T decay) are given in Figs. 
(1) (-'I)- 

Using the parametrization of Eq. (2), for the P — > f £u e differential decay width, we 

get 

2 9 \ 2 



dT 



A 



n 2 it/ |2 \ \I2( 2 2 2\ | Q m £ 

d?-19to*m% G ^ A/ (— /cOf 



x 



+ 



(2q 2 + mj) 



2 

q' 



f + (q 2 ) (2m 2 P + 2m 2 o - q 2 ) + 2(m 2 P - m%) Re[f + (q 2 )r_(q 2 )} 



(q 2 + 2m 2 ) 
q 2 



f + {q 2 ) (m P -m%y 



+ 2(m 2 P - m 2 fo )q 2 Re[f + (q 2 )r.(q 2 )} + f-{q 2 ) 



q 



(19) 



where 



A 



cos 2 for D s — > and , 

for £) -> f £p t and 5^ -> / ^P^ , 
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\V ub \ = 4.31 x 10- 3 for B u -> f £ue , 

Vij = { \V CS \= 0.96 for D s -> , 

1^1 = 0.23 for L> -> / £z7, , 

Taking into account the q 2 dependence of the form factors /+ and /_ and performing 
integration over q 2 and using the lifetimes of B u , D s and D mesons, we get the following 
values for the branching ratios when \ s = 0.19 GeV 2 : 



B(B U ^ 




B(B U ^ 




B(B U - 




B(D S ^ 




B(D S - 




B(D -> 




B[D - 


> foev e ) 



sin 2 e 
2 

sin 2 9 
2 

sin 2 fl 



x (7.17 x 10~ 5 ) , 
x (2.1 x 10" 4 ) , 
x (2.1 x 10~ 4 ) , 



cos 2 # x (3.85 x 10~ 3 
cos 2 # x (4.07 x 10" 3 
sin 2 # 



2 

sin 2 9 



x (4.98 x 10" 4 ) , 
x (5.29 x 10" 4 ) . 



(20) 



We see from (20) that the ratios of the widths 

B(D - fa 



Ri = 
R 2 = 



B(D S 
B(B U 



BID. 



fotoe) ' 



are directly related with the mixing angle 9. On the other hand, as far as the flavor structure 
of /o(980), as is given in Eq. (1), is considered, the ratio 



B(B U -> f £u e 

B(D - f ep e ] 



is independent of the mixing angle 9. Therefore, experimental measurement of the branching 
ratios of B u — > foi^e, D s — > fo£p£ and D — > foive decays can give direct information about 
the mixing angle 9, as well as, about the flavor structure of /o(980) meson. 

In conclusion, we study the semileptonic decay of pseudoscalar mesons to the scalar 
/o(980) meson. The transition form factors are calculated using 3-point QCD sum rule 
analysis and then we estimate the corresponding branching ratios. 
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Figure captions 



Fig. (1) The dependence of the form factor /+ on q 2 at M 2 
m 2 = 0.8 GeV 2 , s = 33 GeV 2 and s' = 1.6 GeV 2 , for the B u 

Fig. (2) The dependence of the form factor /_ on q 2 at M 2 
m 2 = 0.8 GeV 2 , s = 33 GeV 2 and s{, = 1.6 GeV 2 , for the 5 U 

Fig. (3) The dependence of the form factor /+ on q 2 at M 2 
m 2 = 0.6 GeV 2 , s = 7 GeV 2 and s' = 1.6 GeV 2 , for the D s - 

Fig. (4) The dependence of the form factor / + on q 2 at M 2 
m 2 = 0.8 GeV 2 , s = 6 Ge\/ 2 and s' = 1.6 Ge\/ 2 , for the D - 



= 15 GeV 2 , M' 2 = 2 Ge\/ 2 , 
-> / (980)£z7, decay. 

= 15 Gel/ 2 , M' 2 = 2 Gel/ 2 , 
-> / (980)£z7, decay. 

= 6 GeV 2 , M' 2 = 2 Ge\/ 2 , 
■ / (980)^ decay. 

= 6 GeV 2 , M' 2 = 2 GeV 2 , 
/o(980)^ decay. 
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